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In this article we use " Scenario Machine" - the population synthesis simulator- 
to calculate the evolution of the populations of the selected types of X-ray 
sources after the starformation burst with the total mass in binaries 1.5 • lO^M© 
during the first 20 Myr after the burst. Sources of the four types were calculated: 
transient sources- accreting neutron stars with Be- stars; accreting neutron stars 
in pair with supergiants; Cyg X-l-like sources- black holes with supergiants; su- 
peraccreting black holes. We used two values of the a — coef ficient in the 
mass-function: 2.35 (Salpeter's function) and 1.01 ("flat spectrum"). The cal- 
culations were made for two values of the upper limit of the mass-function: 120 
and 30 Mq. For the flat spectrum the number of sources of all types signifi- 
cantly increased. Decreasing of the upper mass limit below the critical mass of a 
black hole formation increase (for the " flat spectrum" ) the number of transient 
sources with neutron stars up to « 300. We give approximating formulae for 
the time dependence of source numbers. 
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1 Introduction 



From the point of view of stellar evolution, objects undergoing bursts of star 
formation are of special interest. One example of such objects are Wolf-Rayet 
galaxies [1]. Considering the only evolution of single stars, Contini et al. [2] 
suggested that some observational propertiesof the galaxy Mrk 712 could be 
understood if the initial mass function were substantially different from the 
standart value. Specifically, they proposed a "flat" mass spectrum with index 
a=l in place of a Salpeter mass spectrum with index a=2.35, with an upper 
mass limit of 12OM0. Later, Schaerer [3], also studying the evolution only of 
single stars, showed that the observations can also be explained using a Salpeter 
initial mass function. 

In [4] , we investigated a burst of star formation under the conditions at the 
center of the Galaxy, and showed that massive X-ray binary systems are some- 
times better indicators of the presence and age of starbursts than the optical 
properties of populations of single stars (color indices, etc.). One of our goals 
here is to draw attentionof population-synthesis specialists to the necessity of 
taking into account the possible effects of synthesis parameters on the evolution 
of populations of binary stars. We consider the evolution of binary systems af- 
ter burst of star formation with two initial mass functions and two upper mass 
limits. 



2 The Model 

Our 'Scenario Machine" population synthesis simulator was first described in 
[5], and recently, a detailed description was given in the review [6] (see also 
[4]). Therefore, here, we will only make note of the parameters used in the 
calculations. 

We considered two power-law initial mass functions, with indices a 2.35 
1.01. 

diss 

^ cx^ (1) 
dM ^ ' 

We calculated 10^ binary systems in each run of the program. The lower 
mass limit for the more massive component was taken to be 0.1 Mq in the 
normalization and IOMq in the computations. Less massive stars do not evolve 
sufficiently over a time < 2 • 10'' years to produce neutron stars and black 
holes (unless their mass is appreciably increased by accretion). The results 
were normalized to a total mass of binary stars in the starburst of 1.5 • 10^ M©. 
We present here calculations for upper mass limits of I2OM0 and 3OM0. We 
chose the rather exotic upper mass limit of 3OM0 for the following reasons. In 
population synthesis modeling, usually two parameters are varied - the slope 
and upper mass limit of the initial mass function (see [11], where, in fact, an 
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upper limit of 30Mq is used). We chose oue smaller value for the upper mass 
limit to be less than the critical mass for the formation of a black hole so that it 
would be possible to study this situation. We took the relative mass distribution 
for the binary components to be flat, i.e., systems with components with similar 
and with very different masses were equally probable. 

We took the major axes of the systems to be in the interval from 10 to lO^i?0, 
with a logarithmically uniform distribution. The magnetic fields of neutron stars 
were assumed not to decay; this is natural, given the relatively short evolution 
period investigated (sec [10] for a discussion of the field decay). A black hole 
was formed when the mass of the pre-supernova star was more than 35 Mq, 70% 
of the initial mass became part of the black hole. The Oppenheimer- Volkov 
lomit was chosen to be 2.5 Mq. We did not use the "enhanced" mass loss 
rate proposed by the Geneva group (see, for example, [3,9]). We will take this 
possibility into account in future papers. 

In anisotropic supernovae explosions, the compact object (neutron star or 
black hole) acquired an additional velocity. These velocities had random di- 
rections and a distribution of amplitudes corresponding to that observed for 
radiopulsars [7], with a characteristic value of 200 km/s (see [8] for a discussion 
of our choice of a lower characteristic velocity than proposed in [7] ) . 

3 Results 

We present our results for four types of sources: (1) X-ray transients associated 
with neutron star - Be star systems; (2) accreting neutron star - supergiant 
systems (one of the types of binary systems that give birth to X-ray pulsars); 
(3) systems with a black hole accreting at a rate above the Eddington limit (SS 
433 may be such an object); and (4) accreting black hole - supergiant systems 
(one example of this type of objects is Cyg X-1). 

The calculated evolution of the number of sources for the two values of the 
initial mass function index are presented in Figs.l (upper mass limit 30 Mq) 
and 2 ((upper mass limit 120 Mq). We have not smoothed the curves, and their 
"noisiness" is due to statistical fluctuations (we calculated 10^ binary systems 
for each mass interval and then normalized the results). X-ray transients were 
the most numerous type of sorce because of the long lifetimes of these systems 
at this stage (the figures show the observable number of systems at each time). 
It stands to reason that we have no information about whether a transient is in 
an active stage at any given moment in time. 

When the upper mass limit is lowered, the number of sources with black 
holes decreases and the number with neutron stars increases. As we decrease 
the upper mass limit below the critical mass for the formation of black holes, 
sources with black holes may not disappear entirely, since they can form if the 
masses of neutron stars are increased to the Oppenheimer - Volkov limit by 
accretion. This is most likely when the accretion is supercritical (and depends 
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on the assumptions about the supercritical accretion stage), as was clear from 
our calculations. For this reason, the number of sources such as SS 433 rose with 
time when the upper mass limit was 30Mq, rather than undergoing the sharp 
decrease observed for the upper mass limit of 120Mq, when the vast majority of 
black holes are formed from massive stars in the first 3-4 million years after the 
starburst. Systems such as Cyg X-1 are not shown for the upper mass limit of 
30Mq, since not a single system of this type formed over the calculation time. 

We present approximation formulas for the convenience in determining the 
number of various types of sources formed in starbursts with arbitrary masses. 
In these expressions, the time t is in million of years. In the case of a Salpeter 
{a = 2.35) initial mass function with an upper limit M„p = 12OM0, we find for 
X-ray transients at times from 5 to 20 million years after the starburst 

N{t) = -0.14 • -I- 5.47 • t - 14.64. (2) 
For superaccreting black holes at times from 4 to 20 million years 

For sources such as Cyg X-1, we have at times from 4 to 20 million years 

= (4) 

We find for binary systems consisting of an accreting neutron star and a 
supergiant at times from 5 to 20 million years 

N{t) = 2.12 • IQ-"' • - 9.6 • 10"'^ • + 0.13 • t - 0.47. (5) 

For a flat initial mass function with an upper mass limit M„p = 12OM0, we 
obtain for X-ray transients at times from 3 to 7 million years 

N{t) = -8.9 • + 1.2 . 10^ . t - 3 • 10^ (6) 
and at times from 7 to 20 million years 

N{t) = -2.8- 1 + 1.2 -10^. (7) 
For superaccreting black holes at times from 4 to 20 millino years, we have 

-.T/x 39.97 

For sources such as Cyg X-1 at times from 4 to 20 million years, 
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Finally, for accreting neutron star - supergiant binaries, we find for times 
from 5 to 20 million years 

N{t) = 1.45 • 10-3 • - 5.96 • 10-2 • + 0.74 • t - 2.41. (10) 

4 Conclusion 

In [4] , we showed that the absolute and relative number of massive baniry sys- 
tems with black holes and neutron stars can serve as a good indicator of the age 
of a burst of star formation. Here, wc considered the evolution of four popula- 
tions of close binary systems with neutron stars or black holes after a starburst 
with a total mass of 1.5 • 10^ in binary systems for the cases of Salpeter 
(a=2.35) and flat (q:=1.01) initial mass functions, following the evolution for 
2 • 10^ years after the starburst. We included the exotic upper mass limit of 
SOMq in our study in order to model populations with a strongly reduced num- 
ber of black holes (note that this situation could also come about in the presence 
of a strong stellar wind associated with a high mass loss rate). 

The transition to the initial mass function index a = 1.01 proposed in [2] 
leads to a signiflcant increase in the number of X-ray sources. Wc expect the 
number of single accreting neutron stars and black holes, and also the number 
of other binary systems, to increase as well, but we have not considered these 
objects here. 

Lowering the upper mass limit below the critical mass for the formation of 
black holes leads to the virtually complete disappearance of sources with black 
holes, but increases the number of transient sources with neutron stars (this 
increase is especially strong for a fiat initial mass function, to « 300 neutron 
star systems). Thus, this change of the initial mass spectrum brings about a 
substantial change in the luminosity of the galaxy at standart and hard X-ray 
energies. In other words, a burst of star formation with a flat initial mass 
function should be accopmanied by a large X-ray flux. The absence of a large 
X-ray luminosity suggests that the initial mass spectrum for the given starburst 
was not flat. 
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